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The question of the exact nature of the phase transition in two-flavor QCD is still under 
discussion. Recent results for small quark masses in simulations with 2 + 1 flavors show scaling 
behavior consistent with the 0(4) or 0(2) universality class. For a precise determination, an 
assessment of deviations from the ideal scaling behavior due to finite quark masses and finite 
simulation volumes is necessary. 

We study the scaling behavior at the chiral phase transition with an effective quark-meson 
model. In our Renormalization Group approach, the quark masses in the model can be varied 
from the chiral limit over a wide range of values, which allows us to estimate scaling deviations 
due to large quark masses and the extent of the scaling region. We conclude that scaling 
deviations are already large at pion masses of 75 MeV, but that the effect is difficult to see in 
the absence of results for even smaller masses. Comparing results only in a narrow window 
of pion masses leads to the observation of apparent scaling behavior. While the scaling devia- 
tions are not necessarily universal, we expect that this may affect current lattice simulation results. 
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By placing the system in a finite box, we investigate the transition between infinite-volume scaling 
behavior and finite-size scaling. We estimate that finite-size scaling behavior can be tested in 
regions where M:^L w 2 — 3, which is smaller than in most current lattice simulations. We expect 
that finite-volume effects are small for pion masses of Mj^ ~ 75 MeV and lattice aspect ratios with 
TL > 8, but that they will become significant when pion masses in lattice simulations become 
smaller. 
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1. Introduction 

The order and the exact nature of the chiral phase transition in QCD are cuiTently being ex- 
plored in lattice simulations [1] which use critical scaling behavior for the analysis [2, 3]. For 
Nf = 2 massless quark flavors, a second-order phase transition in the 0(4) universality class is ex- 
pected, and the more realistic case of A/f = 2 + 1 flavors with two light and one heavy flavor is 
likely also governed by this critical behavior. There have been numerous investigations of critical 
behavior and scaling in the 0(4) universality class [4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14], in spin model 
systems as well as in effective models for chiral symmetry breaking. Most of these studies have 
focussed on exploring the actual critical scaling regime. This requires to effectively go to very 
small masses for the Goldstone modes. 

The idea of our study [15] is to use a simple model for the chiral phase transition to investigate 
the influence of a finite quark mass or a finite volume on the scaling behavior. The model shows 
the same universal behavior which is expected in QCD in the phase transition region. In the model, 
it is easy to obtain results both in the infinite- volume limit in the scaling region and for large pion 
masses and small volumes, where deviations from the scaling behavior become apparent. As far as 
these deviations are mainly determined by the long-range fluctuations, the results are applicable to 
the universal behavior at the chiral transition in QCD. In this way they can shed light on the scaling 
deviations one should expect in the analysis of QCD lattice simulations. 

In the following, we will first present the model that we are using and discuss the scaling 
behavior in the ideal case of small pion masses and lai^ge volumes. We then turn to deviations from 
this scaling behavior in a finite volume. For very small volumes compared to the scale set by the 
pion mass, finite size scaling can be used for the analysis. 

2. The quark-meson model 

The quark-meson model belongs to the family of Nambu-Jona-Lasinio type models and can 
be used to describe the dynamical breaking of the chiral SU (MOl x SU (Nf) r symmetry. We use the 
model with Nf = 2 light flavors, where the chiral symmetry is implemented as an 0(4) symmetry 
in the meson sector. It does not include any gluonic degrees of freedom and is not confining for the 
quarks. We define the model at an ultraviolet scale A by its bare effective action 

rA[^,»F,(/>] = ld^x{^'{i^ + g{a + ir-7ir5))'i'+^{d^(l>f + UA{f)-Ha} (2.1) 

with an effective potential specified by two parameters. In the ground state, the symmetry is dy- 
namically broken, and the first component of = (a, tt) acquires a finite expectation value. There 
are 3 Goldstone bosons associated with the remaining 0(3) symmetry of the ground state. The 
symmetry-breaking parameter H controls the amount of explicit symmetry breaking and the mass 
of the Goldstone bosons. With our choice of parameters, the model undergoes a phase transition at 
Tc ~ 145 MeV to a phase with restored symmetry. 

In order to account for the critical fluctuations at the phase transition, which determine the 
behavior we want to study, we use a non-perturbative Renormalization Group (RG) method and 
solve the Wetterich equation [16] for the RG flow in the local potential approximation. For details 
of our model and calculation, we refer the reader to [15]. 
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Figure 1: Left panel: Longitudinal suszeptibility Xa as a function of temperature f for various pion masses 
Mji < 0.9 MeV. Right panel: Rescaled susceptibility as a function of z for the same pion masses. For these 
small pion masses, all curves fall onto the scaling function with good accuracy. 



Because of the effect of critical fluctuations, close to the critical point the behavior of ther- 
modynamic quantities can be expressed in terms of power laws with universal critical exponents 
which depend on symmetry and dimensionality of the system. As an additional consequence, the 
dependence on temperature T and symmetry-breaking parameter H can be expressed in terms of 
a scaling function of a single scaling variable z- In the case of the order parameter M (the chiral 
condensate or pion decay constant fji in our case), this takes the form 



M{t,h)=h'/^fM{z), z = t/h 



(2.2) 



In order to be able to compare different systems in the same universality class, one defines the 
universal quantities t = {T — Tc)/To and h = H/Hq, where the normalization constants Tq and Hq 
are non-universal and like the critical temperature Tc system-specific. The corresponding relation 
for the longitudinal susceptibility is given by 



/m(z)-JA(z) 



(2.3) 



We have checked carefully that our results for very small pion masses {Mji < 1 MeV) satisfy these 
scaling relations to a very high accuracy. In particular for the susceptibility it is impressive how 
well the results for different pion masses over several orders of magnitude collapse onto the scaling 
function after rescaling, see Fig. 1 and Ref. [15]. 



3. Scaling for large quark masses 

While the expected intinite- volume scaling behavior is realized impressively well in the model 
for a small amount of symmetry breaking, we find that it is necessary to remain close to the critical 
point in order to keep the scaling deviations small. An important question is how large these 
deviations are for realistic values of the pion mass and for pion masses used in lattice simulations. 
Scaling deviations are of course not universal and can reflect the different short-range physics of 
the model compared to QCD. However, the absolute scales (transition temperature) for the effective 
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Figure 2: Left panel: Rescaled d-suszeptibility as a function of z for various pion masses. Right panel: 
Rescaled pion decay constant /^j (order parameter) as a function of z for various pion masses. In the absence 
of the results for small pion masses, the curves for large pion masses appear to fall on the scaling function. 



degrees of freedom at the transition are very similar, and therefore we expect similar deviations, 
provided the behavior is indeed mainly governed by the low-energy degrees of freedom. 

In Fig. 2 we show results for the rescaled susceptibility and the rescaled order parameter as 
a function of the scaling variable z for pion masses of M;i = 75, 138,200 MeV. We find that the 
rescaled results differ significantly from the scaling function that we obtained in the limit of small 
pion masses. The resulting scaling functions are suppressed compared to the proper one, and the 
position of the peak in the susceptibility is shifted compared to the small-pion mass limit. 

In the absence of the results for small pion masses, the rescaled results obtained for large 
pion mass would appear to agree quite well and to actually show scaling behavior: The curves are 
almost on top of one another and differ significantly only far from the critical temperature, below 
the transition. Using these results to determine new values for the normalization constants Tq and 
Hq, they can be rescaled to almost agree with the original scaling function. 

We conclude that the scaling deviations at pion masses of Mji = 75 MeV are already quite 
large, for both the longitudinal susceptibility and the order parameter, but that this effect is difficult 
to observe in the absence of results in the limit of very small symmetry breaking. Comparing only 
results in a narrow window of large pion masses leads to the observation of apparent scaling, while 
these results are still very different from the correct scaling function. 

4. Scaling in a finite volume 

The RG calculation can easily be continued to a finite volume V = [17, 18, 19, 13]. In a 
finite volume, the boundary conditions in the spatial direction for quark fields can be chosen freely; 
we chose periodic boundary conditions, which is the choice taken in most QCD lattice simulations. 
Our main interest in this part of the study was to find out how much the results for the scaling 
behavior change due to the finite volume if we performed the same scaling analysis as in infinite 
volume, i.e determined the scaling function f^iz) for the susceptibility. 

The results are shown in Fig. 3 for a pion mass Mj[ = 75 MeV. In the left panel, we show the 
results for the rescaled susceptibility h^^^^^HoXa as a function of z for several different volume 
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sizes. Significant deviations occur only for volume sizes smaller than L = 8 fm. However, in usual 
lattice simulations at finite temperature, the aspect ratio of the spatial and temporal directions of the 
Euclidean volume is kept fixed and both volume size and temperature are varied at the same time. 
In addition to the results already given in Ref. [15], we show in the right panel of Fig. 3 results 
for this scenario, where the aspect ratio TL of the Euclidean volume has been kept fixed. We find 
that for this pion mass, deviations appear only below TL = 6. We conclude from these results 
that finite-volume effects presumably do not play a large role in current lattice simulations of the 
scaling behavior. However, they are significant and likely to become important for simulations at 
smaller pion masses. 

A somewhat surprising result is the increase in the susceptibility in an intermediate volume 
range that can be seen in Fig. 3, which we attribute partially to quark effects [18]. Naively one 
expects X(y ~ L^, but this behavior is only realized for very small volumes. These finite- volume 
effects are highly dependent on the boundary conditions chosen for the quark fields [18] and can 
also influence the response of the system to a finite chemical potential [20, 21]. 



5. Finite-size scaling 

For very small volume size, when the correlation length of the critical fluctuations is of the 
same order as the volume extent L, scaling with the volume size can be observed. Starting from the 
hypothesis that the system behavior should be unaltered when both system size L and correlation 
length t, are changed such that t,/L remains constant, it is possible to derive scaling functions for 
the volume dependence [4, 13]. These scaling functions will depend on both the infinite- volume 
scahng variable z and a new scaling variable involving L. For example, for the longitudinal suscep- 
tibility one finds 

Xa{t,KL)=Ul''Q^{z.hL^'''') (5.1) 

The new finite-size scaling function Qj^ satisfies the constraint that the infinite- volume scaling 
behavior must be recovered in the infinite-volume limit. 
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Figure 4: Finite-size scaling for the longitudinal susceptibility Xa for z = 0. The unsealed susceptibility 
is shown as a function of the symmetry breaking parameter li (left panel) and the rescaled susceptibility 
is shown as a function of the scaling variable hlP^I^ for various values of the box length L. 



L [fm] 


2 


4 


6 


8 


10 


20 


30 


[MeV] 


308 


139 


85 


60 


45 


19 


11 




3.12 


2.82 


2.59 


2.43 


2.30 


1.94 


1.75 



Table 1: Infinite-volume pion masses M:^{T 0,L o°) which correspond to the value of the scaling 
variable hlP^I^ at the point where the scaling curve in Fig. 4 bends (for z = 0). The dimensionless product 
M-nL contains only long-range quantities and thus provides a benchmark for comparison to QCD simulations. 

This scaling behavior for the longitudinal susceptibility is shown in Fig. 4 for T = Tc. In this 
case the symmetry-breaking parameter h controls the coiTclation length, an increase in h translates 
to larger masses of the fluctuations and smaller con^elation length. In the left panel, it is evident 
that the susceptibility deviates from the asymptotic infinite-volume behavior at different /i-values, 
depending on the volume size. The right panel shows the rescaled susceptibihties, which all fall 
onto the same scaling function. 

Since the volume size is effectively an additional relevant coupling controlling the critical 
behavior, L needs to be renormalized with a non-universal factor Lq, just as T and H, in order to 
facilitate a comparison of different systems. In the present case, it appears more useful to translate 
the scale into terms of the dimensionless quantity M^L. The values of this quantity at the bend, 
where the scaling function for ;^CTbecomes flat, ai^e given in Tab. 1 . In order to observe finite-size 
scaling effects and to use finite-size scaling as a tool for the analysis, it appears necessary to explore 
the region in which MjiL « 2 — 3. This is much smaller than the value achieved in most present 
lattice simulations, where generally Mj^L > 4 — 5 and finite-volume effects are minimized. For a 
systematic exploration of finite-size scaling effects, smaller lattice sizes would be required. 

6. Conclusions 

We have presented results from a scaling analysis of the chiral phase transition in a quark- 
meson model. In such a model the scaling region at small pion masses can be explored as well 
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as deviations from this behavior for large pion masses and finite volumes. In order to properly 
account for critical fluctuations, we used a non-perturbative RG method for the calculation. As 
a main result, we find that for Mj^ > 75 MeV significant deviations appear. This suggests that 
such scaling deviations could also be significant in current lattice QCD simulations. In contrast to 
the effects of finite pion masses, we find finite-volume effects to be small for current simulations 
parameters. In order to use finite- size scaling for the analysis of the phase transition behavior, it is 
necessary to go to much smaller volumes for current pion mass values, such that M^L = 2 — 3. 
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